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0.  Abstract.  Sampling  procedures  are  considered  for  monitoring  the  out- 
put  of  a  production  process.  Sampling  rates  are  allowed  to  depend  on  one  of 
several  levels  of  acceptance  or  rejection,  sampling  less  frequently  when  the 
process  is  in  control.  The  problem  is  formulated  as  a  simple  Markov  process 
whose  properties  yield  the  expected  values  and  the  variances  of  the  number  of 
samples  and  the  numbers  of  various  types  of  adjustments  to  the  process.  Com¬ 
putations  are  given  in  support  of  the  economic  design  of  variable  sampling 
policies  of  this  type. 
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1.  Introduction.  Two  major  facets  of  statistical  quality  control  are 
acceptance  sampling  plans  for  inspecting  lots  of  manufactured  items  and  stat¬ 
istical  control  charts  for  monitoring  the  output  of  a  production  process. 
Central  to  the  present  study  is  the  notion  that  less  intensive  sampling  is 
needed  when  consecutive  lots  of  high  quality  are  produced  or  when  the  produc¬ 
tion  process  continues  to  remain  in  control.  This  idea  is  not  new.  Sampling 
plans  for  the  percent  defective  were  devised  along  these  lines  by  Dodge  (1943) , 
Wald  and  Wolfowitz  (1945),  and  Dodge  and  Torrey  (1951),  who  gave  rules  for 
switching  from  100Z  sampling  to  sampling  a  fraction  f.  Multi-level  inspec¬ 
tion  plans  having  sampling  fractions  in  geometric  progression  were  studied 
by  Lieberman  and  Solomon  (1955)  and  by  Derman,  Littauer  and  Solomon  (1957) 
using  various  rules  for  shifting  to  tightened  inspection  levels.  These  ideas 
carry  over  to  monitoring  production  processes.  In  the  related  problem  of 
monitoring  stream  pollution,  Arnold  (1970)  developed  variable  sampling  poli¬ 
cies  with  subsequent  sampling  rates  depending  on  the  current  level  of  pollu¬ 
tion.  Unlike  process  control,  this  model  makes  no  provision  for  initiating 
corrective  action.  These  policies  were  studied  further  by  Crlgler  (1973)  and 
by  Smeach  and  Jernlgan  (1977).  The  foregoing  procedures,  all  adaptive  in  that 
subsequent  sampling  rates  depend  on  current  outcomes,  effect  greater  economy 
as  the  quality  consistently  improves. 

Data-dependent  sampling  strategies  may  be  described  stochastically. 
Lieberman  and  Solomon  (1955)  thereby  expressed  the  average  fraction  inspected 
in  their  multi-level  plan  in  terms  of  the  steady-state  probabilities  of  a 
Markov  chain.  Similarly,  Arnold  (1970)  found  the  mean  and  Smeach  and  Jernlgan 
(1977)  the  variance  of  the  sample  size  in  Arnold's  variable  sampling  procedure. 

In  this  study  we  develop  adaptive  sampling  policies  for  use  with  control 
charts  for  monitoring  the  output  of  a  production  process.  We  go  beyond  known 


variable  sampling  models,  allowing  for  various  levels  of  acceptance  and  warning 
and  providing  for  various  types  of  rejection.  One  of  several  types  of  adjust¬ 
ments  is  initiated  when  the  process  is  not  in  control.  These  features  enable 
us  to  extend  the  conventional  usage  of  control  charts  to  include  diagnostic 
capabilities  in  multiparameter  cases.  An  example  is  cited.  Following  Arnold 
(1970) ,  we  achieve  variable  sampling  rates  by  varying  the  time  delay  between 
successive  inspections.  On  identifying  our  model  as  an  aperiodic  recurrent 
Harkov  chain,  we  obtain  the  expected  values  of  the  sample  size  and  the  numbers 
of  various  types  of  adjustments  to  production.  Further  studies  yield  exact 
and  approximate  expressions  for  the  variances  of  these  quantities.  Attention 
then  is  given  to  the  design  of  time-delay  sampling  policies. 

2.  The  Model  and  Assumptions.  Let  Y(t*l)  be  the  typical  quality  char¬ 
acteristics  of  a  production  process,  taking  values  in  the  t-dimensional  Eucli¬ 
dean  space  R1"  according  to  a  probability  measure  P(*).  We  partition  RC  into 
regions  (1^,  ...  ,1^,  R^,  ...  ,Rm}  consisting  of  outcomes  of  two  types.  Out¬ 
comes  in  {1^,  ...  ,1^1  indicate  different  levels  of  acceptance  and  warning; 
the  outcomes  in  {R^,  ...  ,Rm)  indicate  different  types  of  rejections  requiring 
different  types  of  adjustments  to  production  when  the  process  is  not  in  control. 

To  fix  ideas  for  the  case  k  *  2  and  m  =  1,  consider  the  use  of  Hotelling's 
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T  chart  (cf.  Hotelling  (1947))  for  monitoring  the  means  of  a  bicharacteristic 
process  when  the  process  dispersion  parameters  are  known.  The  regions  {1^, 

I^,  R)  may  be  constructed  as  in  Figure  1  with  1^  an  acceptance  region,  a 
warning  region,  and  R  a  rejection  region  signalling  that  the  process  is  not 
in  control  and  that  corrective  action  is  needed.  An  example  having  k  *  1  and 
m  *  3  arises  in  monitoring  the  means  of  a  bicharacteristic  process  against  one¬ 
sided  upper  alternatives,  the  object  being  to  diagnose  which  variables  require 


adjustment  when  the  process  is  not  in  control.  The  regions  (I,  R^,  } 

as  in  Figure  2  may  be  chosen  with  I  an  acceptance  region,  R^  a  rejection  region 
requiring  adjustment  to  but  not  u2,  R2  a  re8ion  requiring  adjustment  to  p2 
but  not  and  R^  a  region  requiring  adjustments  to  both  p^  and  p2* 

We  henceforth  assume  that  items  yielding  the  characteristics  Y  are  pro¬ 
duced  in  sequence,  one  unit  of  time  having  elapsed  in  the  production  of  each. 
Let  {d^,  ...  ,d^}  be  positive  delay  parameters  not  exceeding  some  maximum  delay 
A,  and  suppose  for  each  j  =  1,  2,  ...  ,m  that  adjustments  originating  in  the 
region  R^  require  A^  time  units  for  completion.  Our  basic  sampling  plan  is  the 
following. 

Time-Delay  Inspection  Plan.  If  Y  t  1^,  delay  d^  time  units  before  in¬ 
specting  the  next  item.  If  Y  e  Rj »  begin  adjustment  as  appropriate,  delay 
sampling  for  Aj  time  units,  and  then  resume  inspection  of  the  adjusted  process. 

Note  that  the  basic  sampling  plan  is  advantageous  when  inspection  costs 
are  high  and  it  is  more  economical  to  forego  further  inspection  during  adjust¬ 
ment  and  to  sell  at  reduced  prices  the  items  of  questionable  quality.  A  modi¬ 
fication  of  the  basic  plan  is  to  inspect  all  items  consecutively  for  A^  -  1 
time  intervals  during  adjustment,  and  then  to  resume  inspection  of  the  adjusted 
process.  This  modification  may  be  appropriate  when  inspections  are  critical, 
ensuring  as  it  does  that  all  items  of  questionable  quality  will  have  been  in¬ 
spected. 

In  order  to  study  the  stochastic  behavior  of  our  sampling  policies  effi¬ 
ciently,  we  approach  the  problem  as  an  application  of  the  theory  of  discrete 
Markov  chains.  The  latter  are  characterized  in  terms  of  their  possible  states 
and  the  matrix  of  one-step  transition  probabilities.  This  program  is  carried 
out  in  the  following  sections,  first  for  a  single  rejection  model,  then  for 


FIGURE  2.  Control  limits  for  monitoring  the  means 
of  a  bicharacteristic  process  against  one-sided 
alternatives  with  diagnostic  capabilities. 
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the  general  case.  The  findings  are  applied  to  the  design  of  variable  sampling 
plans  of  the  types  considered  here. 

3.  Single  Rejection  Models.  In  sequence  we  identify  the  Markovian  states; 
we  give  the  one-step  transition  probabilities  in  terms  of  the  natural  parameters 
of  the  problem;  and  we  develop  some  of  the  stochastic  properties  of  our  time- 
delay  monitoring  scheme.  Proofs  are  given  in  detail  in  preparation  for  the 
sketches  outlined  in  the  section  following. 

Consider  the  regions  {1^,  •••  ,1^*  R)  with  respective  probabilities  {p^, 

...  ,p^  q)  under  P(*)  such  that  p^  +  . . .  +  p^  +  q  «  1.  An  outcome  is  said  to 
be  acceptable  if  Y  e  1^  for  some  i,  and  to  be  unacceptable  if  Y  e  R,  in  which 
case  adjustments  requiring  A  time  units  are  initiated.  Let  be  the  sampling 
state  that  the  most  recent  sample  was  taken  i  time  units  previously  and  Y  was 
found  acceptable,  where  1  <  i  <  6  and  6  =  maxCd^,  ...  .d^},  and  let  t  denote 
the  event  of  being  in  sampling  state  E^^  at  some  nonnegative  time  t.  Similarly 
let  Fj  be  the  sampling  state  that  Y  was  found  unacceptable  and  adjustment 
commenced  j  time  units  previously,  where  1  <  j  <  A.  Clearly  {E. ,  ...  ,E_, 

F^,  ...  ,F^}  are  the  exclusive  and  exhaustive  states  of  our  time-delay  sampling 
scheme.  The  procedure  is  initiated  by  taking  a  sample  at  time  t  =  0  when  the 
process  is  in  control,  the  sampling  procedure  being  in  state  E^  at  time  t  =  1. 

It  is  important  to  note  that  the  sampling  state  occupied  at  time  unit  t  is  de¬ 
termined  before  any  sample  is  taken  and  before  any  adjustments  to  the  process 
have  begun. 


The  matrix  of  one-step  transition  probabilities  is  found  in  terms  of  the 

parameters  of  the  problem  as  follows.  Retain  the  order  fE, ,  ...  ,E. ,  F, ,  ...  , 

1  ol 

F.l,  and  consider  the  corresponding  stochastic  matrix  P  =  [P  ]  of  order 
a  v  rs 

(<5+A)  x  (6+A).  For  1  <  t  <  5  let  ir  be  the  probability  of  delaying  exactly 
t  time  units  before  the  next  inspection  given  that  the  most  recent  inspection 
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is  acceptable.  Because  inspection  is  delayed  exactly  d^  time  units  when  Y  e  1^ 
and  this  occurs  with  probability  p^,  it  follows  that 


it 


t 


p^  if  t  =  d^  for  some  i  =  1,  2,  ...  ,k 
0  otherwise. 


(3.1) 


In  one  time  unit  the  only  possible  state  transitions  are  (i)  E..  to  Ej+^,  i  =  1, 

2,  ...  ,6-1;  (ii)  E^  to  E^,  i  =  1,  2,  ...  ,6;  (iii)  E^  to  F^,  i  =  1,  2,  ...  ,6; 

(iv)  to  Fj+1»  j  =  1,  2,  ...  ,A-1;  (v)  F^  to  E^;  and  (vi)  FA  to  F^.  The  only 

nonzero  elements  of  P  accordingly  are  those  corresponding  to  these  six  state 

transition  types.  The  structure  of  P  is  reflected  in  Table  1;  its  typical 

elements  as  functions  of  {it.,  ...  ,ir.,  q}  are  considered  next  on  a  case-by-case 

X  o 

basis.  To  illustrate  the  types  of  arguments  used  note  that 


Pi,i+1."  P(Ei+l,t+llEi,t) 


=  1  -  P (inspect  at  the  next  time  unit|E.  ) 

x  *  t 

^  P (inspection  delays  exactly  i  time  units) 

P (inspection  delays  at  least  i  time  units) 


Now  let  ir(r)  =  (tt^  +  irr+^  +  •••  +  it^)  be  the  probability  that  inspection  is 
delayed  at  least  r  time  units  following  a  favorable  inspection.  Routine  argu¬ 
ments  are  used  to  establish  the  following. 

Case  (i).  P^  i+1  =  1  -  t^/ttU)  ,  i  =  1,  2,  ...  ,6-1; 

Case  (ii).  P^  ^  =  (l-qH^/ird) ,  i  =  1,  2,  ...  ,6; 

Case  (iii).  P.^  =  qir±/ir(i),  i  *  1,  2,  ...  ,6; 

Case  (iv) .  6+i+l  =  ^  ~  2,  ...  , A— 1; 

Case  (v).  P6+A>1  =  1-q; 

Case  (vi).  Pm>6+1  =  q; 

Case  (vii) .  P^j  =  0,  otherwise. 


TABLE  1.  The  one-step  transition  probabilities  Prg  from  state  r  to  state  s 
in  a  time-delay  sampling  scheme  with  single  rejection  region. 
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It  is  seen  from  the  one-step  transition  matrix  P  that  the  Markov  process 
is  aperiodic  and  recurrent.  Its  properties  yield  insight  into  the  stochastic 
behavior  of  the  number  of  inspections  and  the  number  of  adjustments  to  the  pro¬ 
cess  in  monitoring  over  a  period  of  length  T.  It  is  readily  seen  that  the  num¬ 
ber  L^,  of  acceptable  inspections  in  a  time  interval  (0,T]  is  equal  to  the  number 
of  visits  to  state  in  that  interval.  Similarly  the  number  MT  of  unacceptable 
inspections  is  equal  to  the  number  of  visits  to  state  in  the  interval  (0,T], 
and  this  is  precisely  the  number  of  adjustments.  Let  P^”^  be  the  (i,j)  element 
of  the  nth  power  of  P.  Then  is  the  probability  of  a  favorable  inspection 

at  time  t,  i.e.,  the  probability  of  transition  from  to  E^  in  t  steps.  Sim¬ 
ilarly  P^g+i  is  the  probability  that  adjustment  is  initiated  at  time  unit  t, 
i.e.,  the  probability  of  transition  from  E^  to  in  t  steps.  It  follows  using 

standard  arguments  (cf.  Arnold  (1970),  for  example)  that  the  expectations  are 

T 


E(U.)  =  £  P(X) 

t=l 


(3.2) 


E«r>  '  ^  • 


(3.3) 


Let  Nt  be  the  total  number  of  inspections  carried  out  in  the  interval  (0,TJ. 
The  mean  E(N^,)  and  the  variances  V(LT)  and  V(NT)  are  given  in  the  following. 


i  * 


THEOREM  3.1.  Let  LT  and  respectively  be  the  number  of  acceptable  and  the 
total  number  of  inspections  in  monitoring  production  over  the  interval  (0,T] 
using  time-delay  sampling  policies.  Then  the  mean  and  variance  of  are  given 


respectively  by 


E(v  ■  j/u  + 


V(Nl>  =  i-1  j-l  a* I J I  '  3iaj] 


(3.4) 

(3.5) 
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where  afc  -  and  a*  -  (l-q)at  +  q(P^fl  +  p£i>6+1>-  The  mean 

and  variance  of  are  given  respectively  by  (3.2)  and 


T  T 


v<v '  i«i  A  bu-ji  -  w 


(3.6) 


where  bt  *  P^^  • 


Proof.  Expression  (3.4)  follows  directly  from  (3.2)  and  (3.3)  and  the  relation 

+  ^T*  Expression  (3.6)  follows  in  the  manner  of  Smeach  and  Jernigan 

(1977)  starting  with  bt,  the  probability  of  an  acceptable  inspection.  To  verify 

(3.5),  note  that  a.  =  pf^  +  is  the  probability  that  an  inspection  is  made 

C  11  1,0+1 

at  time  t.  Following  Smeach  and  Jernigan  (1977),  let  Z^  be  the  indicator  of  this 
event  and  note  that  NT  *  +  ...  +  Z^,.  Using  the  fact  that  the  covariance  of 

Z^  and  Zj  may  be  expressed  as  CovCZ^Zj)  =  P(Zj=l,  Z^=l)  ~  E(Z^)  E(Zj),  we  have 


T  T 

V(N_)  -  Z  Z  (P(Z  -1,  Z  =1)  -  E(Z . )  E(Z  )]  . 
i-1  j-1  1  3  J 


(3.7) 


Let  r  ■  min(i,j)  and  s  ■  |i-j|  and  consider  the  conditional  probability  P(Zr+s*l| 
Zr»l) .  Given  that  Z^  *  1,  an  inspection  was  made  at  time  t  =  r  and  thus  the 
process  is  either  in  sampling  state  Ej  at  t  *  r  +  1  or  In  sampling  state  F^  at 
t  *  r  +  1,  these  events  occurring  with  probabilities  (1-q)  and  q,  respectively. 
The  law  of  total  probability  asserts  that 


W^r"1*  "  P(insPect  at  t  -  r  +  s|El  r+1)  P(E1>r+1) 

+  P( inspect  at  t  ■  r  +  slF^r+l))  P(Fj(r+l)) 

<P11>  +  Pl,6+1)  <1-q)  +  <P6+1,1  +  P6+l,6-»-l)q 


(3.8) 


where  F^(t)  la  Che  event  that  the  process  is  in  sampling  state  F^  at  time  t. 
On  writing 


p(z1-i,  Zj-i)  -  p(zr+s-ijzr-i)  P(Zr-l) 


(3.9) 
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and  noting  that  P(Z  “1)  =  E(Z  )  =  P^  +  ,  we  combine  (3.9)  and  (3.8) 

r  r  xx  x  9 6+x 

with  (3.7)  to  get  (3.5). 

Suppose  the  inspection  plan  is  modified  so  that,  for  any  adjustment  ini¬ 
tiated  in  (0,T],  the  next  A  items  are  inspected  during  adjustment.  Let  N*  be 
the  total  number  of  inspections  in  (0,T]  under  this  scheme.  Then  (3. A)  is  modi¬ 
fied  to 

T  T 

E(N*)  =  I  P.(5}  +  A  £  pf0...  .  (3.10) 

T  t=l  11  t-1  1,6+1 

In  a  manner  similar  to  our  treatment  of  Lp  and  Np,  we  may  investigate  further 
stochastic  properties  of  the  number  Mp  of  adjustments  to  the  process.  For  ease  of 
reference  our  main  results  are  summarized  in  the  following. 


THEOREM  3.2.  Let  Mp  be  the  number  of  adjustments  in  monitoring  production  over 
the  interval  (0,T]  using  time-delay  sampling  policies.  Then  E(Mp)  and  V(Mp)  are 
given  respectively  by  (3.3)  and 


T  T 
V(M  )  =  I  Z 
i=l  j-1 


^Cmin(i,j)  Cfi-j| 


Vi1 


(3.11) 


where  c 


P1,<S+1  and  CT  =  P  5+1, 6+1’ 


Proof.  Note  that  ct  =  is  the  probability  that  an  adjustment  is  initiated 

at  time  t.  Let  be  the  indicator  for  this  event  and  note  that  Mp  *  Zp  +  . . .  + 

Zp.  The  proof  now  parallels  that  of  Theorem  3.1  on  interpreting  the  conditional 

probability  P(Zr+a“l |Zr*l)  suitably.  Given  that  Zr»l,  the  process  was  adjusted  at 

time  t  ■  r  and  thus  is  in  sampling  state  F.  at  t  «  r  +  1.  In  order  that  Z^^  *  1, 

i  r+8 


adjustment  must  be  initiated  at  t  *  r  +  s,  assuring  that  the  process  is  in  sampling 

state  F.  at  t  *  r  +  s  +  1.  Given  that  Z  *  1,  the  event  Z  .  ■  1  thus  is  attained 

i  r  r+s 

through  the  transition  from  sampling  state  Fp  to  state  Fp  in  exactly  s  steps. 
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i.e.,  P(Zr+s-l|Zr=l)  =  a+1.  Because  P(Zr=l)  =  E(Zr>  =  the  expres¬ 

sion  corresponding  to  (3.9)  is 


PCZj-l,  Zj-1) 


[P 


(s) 

6+1, 6+1 


1  [P 


(r) 

1,6+1 


1  • 


(3.12) 


The  remainder  of  the  proof  is  identical  to  that  of  Theorem  3.1. 


The  stochastic  properties  of  the  number  of  inspections  and  the  number  of 
adjustments,  as  set  forth  in  Theorems  3.1  and  3.2,  assume  a  vital  role  in  com¬ 
paring  alternative  time-delay  sampling  plans.  We  return  to  this  topic  in  Sec¬ 
tion  5. 


4.  Multiple  Rejection  Models.  Motivated  in  part  by  the  need  for  diagnostic 
capabilities  mentioned  with  reference  to  Figure  2,  we  consider  the  case  of  mul¬ 
tiple  rejection  regions  requiring  different  types  of  adjustments  to  the  process. 

An  extension  of  the  arguments  of  the  preceding  section  yields  stochastic  proper¬ 
ties  of  the  number  of  samples  and  the  numbers  of  the  several  types  of  adjustments 
in  the  time  interval  (0,T]. 

Consider  the  regions  (l^,  ...  ,1^,  R^,  ...  .R^)  with  respective  probabili¬ 
ties  (p^,  ...  ,p^,  q^,  ...  .q^)  under  P(’)  such  that  P^  +  •••  +  p^  +  q^  +  ...  + 

*  1.  An  outcome  is  said  to  be  acceptable  if  Y  e  1^  for  some  i  ■  1,  2,  ...  ,k, 
and  to  be  unacceptable  of  type  j  if  Y  e  Rj ,  in  which  case  the  process  is  said 
to  be  in  rej ectlon  mode  for  A^  time  units  while  the  process  is  undergoing  adjust¬ 
ment.  Define  the  sampling  states  {E^,  ...  ,Efi}  as  before,  and  let  be  the 
sampling  state  that  the  process  is  currently  in  rejection  mode  i,  having  entered 

that  mode  j  time  units  previously,  for  j  ■  1,  2,  ...  ,A^  and  i  «  1,  2 . m. 

Consider  the  states  (E^,  ...  ,Eg,  F^,  ...  ,F^  ,  ...  .F^,  ...  in  that 

order;  these  clearly  are  the  exclusive  and  exhaustive  states  of  our  time-delay 
sampling  policies  with  multiple  rejection  regions.  Let  P  “  tP^j]  be  the 
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corresponding  matrix  of  one-step  transition  probabilities;  its  structure  and  typ¬ 
ical  elements  are  found  as  in  Section  3. 

In  one  time  unit  there  are  six  possible  state  transitions,  namely,  (i) 
to  E1+1,  i  ■  1,  2,  ...  ,6-1;  (ii)  Ei  to  i  =  1,  2,  ...  ,6;  (iii)  EA  to  F^, 
i  *  1,  2,  ...  ,6,  j  *  1,  2,  ...  ,m;  (iv)  to  Ft  j+1»  i  ■  1,  2,  ...  ,m,  j  » 

1,  2,  ...  ,A^-1;  (v)  FiA  to  E^,  i  =  1,  2,  ...  ,m;  and  (vi)  F^  to  F^,  i  * 

1,  2,  ...  ,m,  j  =*  1,  2,  ...  ,m.  Thus  the  only  nonzero  elements  of  P  correspond 


to  transitions  of  these  types. 

To  evaluate  typical  elements  of  P  let  R  =  u  . . .  u  and  q  *  +  . . .  + 

q^,  and  define  {A(r)  =  AQ  +  A^  +  ...  +  Ar;  r  =  0,  1,  ...  ,m}  with  Aq  ■  0.  Re¬ 
calling  that  if(r)  =»  (it  +  ir  .,+...+  irr) ,  we  proceed  as  in  Section  3  to  eval- 

r  r*ri  o 

uate  Py  in  the  following  typical  cases. 


Case  (i). 


i ,  i+1 


1  -  iTj/ird),  i  -  1,  2,  ...  ,6-1; 


Case  (ii).  ?i  ^  =  (1-q)  ir^/ir(i),  i  *  1,  2 . 6; 

Case  (iii).  P^  ,5+A(j— i)+i  *  if^/’n(i)»  1  *  1»  2,  ...  ,6,  j  *  1>  2,  ...  ,m; 

Case  (iv) .  P 


6+A(i-l)+j,6+A(i-l)+j+l 
Case  (v).  P5+A(1)>1  =  1-q,  i  =  1,  2,  . 


1,  i  =  1,  2, . . .  ,m,  j  -  1,  2, . . .  ,Aj-l; 


.m; 


Case  (vi), 


6+A(i),6+A(j-l)+l 


q j ,  R  =  1,  2,  ...  ,m,  j  *  1,  2, 


Case  (vii).  P^  *  0,  otherwise. 


Let  Lt  be  the  number  of  acceptable  inspections  in  the  time  interval  (0,T] 
and  let  (M^T»  M2T*  ,MmT^  resPectively  be  the  number  of  adjustments  of  types 

1,  2,  ...  ,m.  Clearly  the  number  of  adjustments  of  type  j  initiated  during  the 
period  is  equal  to  the  number  of  entries  into  rejection  mode  j,  and  this  is  pre¬ 
cisely  the  number  of  visits  of  the  process  to  sampling  state  F^  in  (0,T].  Be¬ 
cause  the  process  is  in  sampling  state  Ej  at  time  t  ■  1,  it  follows  as  before 
that  the  probability  of  a  favorable  inspection  at  time  t  is  P^^ .  Similarly 
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i 

i 


S*A(j-l)+l  *8  the  Probabillty  that  Y  Is  unacceptable  of  type  j  at  time  t  and 
thus  that  an  adjustment  of  type  j  Is  initiated.  Expressions  (3.2)  and  (3.6)  for 
E(L,j,)  and  V(Lj.)  continue  to  hold  with  bfc  *  P^jj^  defined  in  terms  of  the  current 
version  of  P^n\  Using  familiar  arguments  we  have 

T 

E<MjT)  "  Pl,$fA(j-l)+l  5  J  “  2*  •••  (4,1) 

Low  order  moments  of  the  distribution  of  the  total  number  of  inspections 
are  given  in  the  following. 


THEOREM  4.1.  Let  be  the  total  number  of  inspections  in  monitoring  production 
over  the  interval  (0,T]  using  time-delay  sampling  with  multiple  rejection  capa¬ 
bilities.  Then  the  mean  and  variance  of  are  given  respectively  by 


e(nt) 


T  (t)  ® 

I  P,,  +  £ 


£  P 


(t) 


t-l  11  j-l  t-l 


(4.2) 


V  (N_) 


'T> '  )h  “Ti-j  i '  aiai' 


(4.3) 


where  at  -  pJ^  +  ^  $+*(•)_!)+!  and 


a*  -  (l-q)at  +  ^  +  ^  P6+A(j-1)+1,6+A(k-1)+11  * 


(4.4) 


f 

K't 


•  t 


Proof .  Expression  (4.2)  follows  directly  from  (3.2)  and  (4.1)  in  view  of  the 
relation  N^,  *  Uj,  +  +  ...  +  M^.  To  validate  expression  (4.3),  let  Z£  be 

an  indicator  for  the  event  that  an  Inspection  is  made  at  time  t  and  note  that 
E(Z^)  -  a£ .  Our  proof  parallels  that  of  Theorem  3.1,  requiring  the  conditional 
probability  P(Zr+#»l|Zr»l).  Given  that  Zr«l,  the  event  Z^ ( ^-1  requires  transi¬ 
tion  from  one  of  the  states  (E^,  Fj^,  •••  tFa^>  at  time  t  ■  r  +  1  to  one  of  those 
same  states  at  t  ■  r  +  i  +  1.  An  extension  of  our  earlier  arguments  yields 


(4.5) 


P(Zr+s-l|Zr=l)  -  P (inspect  at  t-r+s  >r+1)p(E1>r+1) 
m 

+  E  P (inspect  at  t»r+s|F  .(r+l))P(F  (r+1)) 

j-1  2  J 

where  F^(t)  is  the  event  that  the  process  is  in  sampling  state  F^  at  time  t. 

On  evaluating  (4.5)  we  get  (4.4).  The  remainder  of  the  proof  is  identical  to 
the  proofs  of  Theorems  3.1  and  3.2. 

The  remainder  of  this  section  is  concerned  with  the  number  of  adjustments 
of  type  j  undertaken  while  production  is  monitored  over  the  interval  (0,T].  For 
ease  of  reference  the  results  are  summarized  as  follows. 


THEOREM  4.2.  For  j  -  1,  2,  ...  ,m  let  MjT  be  the  number  of  adjustments  of  type 
j  in  monitoring  production  over  the  interval  (0,T]  using  time-delay  sampling 
with  multiple  rejection  capabilities.  Then  the  mean  and  variance  of  are 
given  respectively  by  (4.1)  and 


V(Mjt) 


T  T 
X  l 
u-1  v=l 


^min(u,v)  c  u-vl 


-  bjbj] 


U  V 


(4.6) 


where  bj  -  Ix>$+A(;j_1)+1  and  *  Pfi+A(j-1)+1, 5+A(j-l)+l' 

Proof .  The  proof  parallels  that  of  Theorem  3.2  with  Zfc  an  indicator  for  the  event 
that  a  type  j  adjustment  begins  at  time  t.  The  essential  difference  lies  in  eval¬ 
uating  the  conditional  probability  P(Zj+s-l|Zr-l) .  However,  given  that  Zf-1, 

the  process  is  in  sampling  state  F.,  at  time  t  -  r  +  1,  and  the  event  Z...  -1 

j  i  r+s 

implies  the  sampling  state  F^att-r+s+1.  The  probability  of  transition 
from  state  Fj^  to  state  Fjl  in  exactly  s  steps  thus  gives 


p(zr+8-1! 


Zr-1)  -  p 


(8) 

6+A(j-l)+l,6+A(j-l)+l  * 


(4.7) 


The  remainder  of  the  proof  is  identical  to  that  of  Theorem  3.2. 
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5.  Choice  of  Sampling  Plans.  The  developments  of  the  foregoing  sections 
bear  heavily  on  the  choice  of  adaptive  sampling  plans  for  monitoring  production. 

We  assume  that  the  types  and  number  of  regions  may  be  chosen  in  the  context  of 
eaeh  particular  problem,  but  we  caution  that  the  remaining  parameters  may  be 
excessive  in  number.  These  Include  the  probabilities  {p^,  ...  ,p^,  q^,  ...  , q^}, 
determined  through  choice  of  the  boundaries  of  the  regions,  and  the  delay  para¬ 
meters  {d^,  ...  ,d^}.  A  detailed  treatment  of  economic  models  and  optimal  sam¬ 
pling  plans  is  beyond  the  scope  of  this  study;  various  search  algorithms  may  be 
put  to  this  purpose.  Instead  we  consider  more  elementary  aspects  of  the  choice 
of  sampling  plans,  illustrating  these  numerically  in  selected  cases. 

Arnold  (1970)  and  Crigler  (1973)  focused  mainly  on  the  expected  number  of 
inspections  in  choosing  among  the  time-delay  sampling  policies  considered  there. 
However,  Smeach  and  Jernigan  (1977)  emphasized  also  the  importance  of  the  vari¬ 
ance  of  the  number  of  Inspections  in  making  that  choice.  In  the  model  under 
current  study,  where  adjustments  to  production  may  be  considerably  more  expensive 
than  routine  Inspections,  it  appears  desirable  to  consider  the  mean  and  variance 
of  the  number  of  adjustments  as  well. 

To  fix  ideas  consider  again  monitoring  a  bicharacteristic  production  process 
2 

using  Hotelling's  T  chart  with  acceptance  (Ij),  warning  (Ig),  and  rejection  (R) 
regions  as  in  Figure  1  such  that  p^  =  0,65,  *  0*30,  and  q  ■  0.05.  We  retain 

these  values  and  vary  the  design  through  choice  of  the  delay  parameters  {d^,  d^} 
subject  to  d^  >  d,  and  a  maximum  allowable  delay  of  A  *  10  time  units.  It  is 
assumed  that  adjustment  of  the  process  requires  A  *  5  time  units,  and  that  the 
process  is  to  be  monitored  for  T  ■  1000  time  units.  Eight  different  sampling 
plans  are  listed  in  Table  2  according  to  {d^,  cU } ,  together  with  the  values  of 
E(H^)  and  E(Hj,)  determined  as  in  Section  3.  Suppose  budgetary  allocations  will 
support  a  sampling  plan  averaging  between  155  and  160  inspections.  The  candidates 


TABLE  2.  Expected  numbers  of  Inspections  (NT)  and  adjustments  (Kj.)  under  alter¬ 
native  plans  In  a  time-delay  sampling  model  having  regions  {1^,  R)  with 
probabilities  {0.65,  0.30,  0.05},  operating  over  the  time  interval  (0,  1000] 
with  A  *  5. 


Plan 

dl 

d2 

e[nt] 

EIM,.] 

I 

2 

1 

540.45 

27.02 

II 

9 

1 

155.97 

7.80 

III 

3 

2 

356.84 

17.84 

IV 

8 

3 

157.10 

7.86 

V 

8 

4 

.  149.97 

7.50 

VI 

7 

5 

158.30 

7.92 

VII 

9 

6 

126.14 

6.31 

VIII 


10 


9 


105.35 


5.27 
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are  Plans  II,  IV,  and  VI,  and  these  accordingly  are  examined  In  greater  detail 
In  Table  3  with  regard  to  the  means  and  variances  of  both  NT,  the  number  of 
Inspections,  and  Kj,,  the  number  of  adjustments  to  production.  Computations  were 
carried  out  as  prescribed  in  Section  3.  As  the  means  and  variances  of  Hj,  are 
comparable,  we  distinguish  among  the  three  plans  on  the  basis  of  N^,  noting  that 
its  variance  is  considerably  smaller  for  Plan  VI.  On  these  grounds  Plan  VI 
accordingly  would  be  recommended. 

We  conclude  this  section  with  notes  on  computations.  For  T  as  large  as  1000 
the  computations  of  means  and  variances  are  formidable.  In  addition  to  estab¬ 
lishing  exact  variance  formulas  as  noted  earlier,  Smeach  and  Jernigan  (1977)  gave 
approximations  to  the  mean  and  variance  of  in  the  model  considered  there. 

Their  approximations  carry  over  to  the  present  study  in  the  following  manner. 

Because  the  Markovian  sampling  scheme  of  Section  4  is  aperiodic  and  recur¬ 
rent,  the  limiting  probabilities  for  the  several  states  are  given  by 

^  1  ,-1 


(«*>) 

P1 


[i  +  z  (  n  p4 


i-l  j-1 


j  ,  j+r 


(5.1) 


n-l 


(“>)  (  T>  V  (■>)  » 

"«  "i  •  2>  . * 


(5.2) 


where  y  ■  6  +  A(m).  Then  the  limiting  probability  a  of  making  an  inspection 

oo 

is  given  by 

/_\  m  t  \ 

(5.3) 


(«) 
'1 


a  -  P,'  '  +  E  p 


j-1 


(o») 

6+A(j-l)+l  • 


If  the  convergence  of  a^  to  a^,  is  rapid,  then  instead  of  (4.2)  one  may  use  the 


approximation 


E(N_)  -  Ta 

X  <x 


(5.4) 


To  approximate  variances,  note  from  the  definition  (4.4)  of  a*  that  its 


limiting  value  is 


d 


TABLE  3.  Means  and  variances  of  the  number  of  Inspections  (NT>  and  the  number 
of  adjustments  (M^)  for  three  plans  listed  in  Table  2. 
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a* 


(1-qMpJ*0 


+  Z 

J-l 


P 


(“) 

6+A(j-l)+l 


] 


m  r  (») 

+  £  qJpi 

j=i  3 


+ 


m 

2  P 
k=l 


(“) 

6+A(k-l)+l 


] 


(5.5) 


i.e.,  a*  ■  a^.  For  some  specified  error  e  >  0  let  =  min{t||at  -  aw|  <  e), 

“  min{tj|a*  -  a^J  <  e),  and  t  ■  max{t^,  T2}.  Following  Smeach  and  Jernigan 
(1977),  we  compute 

V«W  '  V(V  ■ 2  A(vWi,  -  Vt+i> ♦  at+i  -  1  • 

i”l  1  1 

Because  a  . ,  -*■  a  and  a*  a  ,  for  t  >  i  we  have  the  approximation 
t+1  00  u  ® 

T 

V(Nt+l)  “  V(V  *  aJl-(2t+l)  aro  +  2  Z  a*]  .  (5.7) 

i=l 


It  follows  that 


T 

V(N_)  A  Ta  (l-(2T+l)a  +2  Z  a*]  .  (5.8) 

T  *  i«i  1 

To  compare  the  approximation  with  the  exact  computation,  we  return  to  the 
three  sampling  plans  given  in  Table  3.  Rather  than  fixing  e  arbitrarily,  we 
specify  t  ■  100  and  apply  the  approximation  (5.8)  with  T  =  1000.  The  results 
are  reported  in  Table  4,  suggesting  that  considerably  more  than  100  terms  would 
be  required  for  these  cases  in  order  to  achieve  acceptable  accuracy. 


6.  Concluding  Remarks.  Several  variations  of  our  basic  model  can  be  given. 

If  no  rejection  is  specified  on  the  range  of  the  quality  characteristics  Y,  we 
obtain  the  sampling  structure  considered  in  Arnold  (1970) ,  Crigler  (1973) ,  and 
Smeach  and  Jernigan  (1977).  The  sampling  states  are  {E^,  ...  ,E^};  nonzero 
elements  of  the  matrix  P(6*6)  of  one-step  transition  probabilities  are  {P^ 
i  ■  1,  2,  ...  ,6}  and  (P.  ... ;  i  -  1,  2,  ...  ,6-1)  with  P.  ,-l;  and  the  expected 

1 9  lTl  0  »  A 

sample  size  in  the  interval  (0,T]  is 

i 


TABLE  4.  Exact  and  approximate  variance  of  the  number  of  inspections  using 
x  =  100  terms  for  three  sampling  plans  listed  in  Table  2. 


Plan  d^  V[N^,]  Approximation 


II 

9 

1 

104.83 

144.09 

IV 

8 

3 

71.39 

109.44 

VI 


7 


5 


55.25 


92.53 
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E(N  )  =  E  P(*} 
t-1  X± 


(6.1) 


Another  Interesting  possibility  is  to  terminate  sampling  when  Y  e  Rj  for  any 
j  *  1,  2,  ...  ,m.  The  sampling  states  are  {E^,  ...  ,E^,  F^,  ...  .F^}  where  F^ 
is  the  absorbing  state  corresponding  to  Y  e  .  Nonzero  elements  of  the  matrix 
of  one-step  transition  probabilities  are  {P^  i  =  1,  2,  ...  ,6},  {P^ 

i  ”  1*  2*  *•'  »6-1^  *Pi,6+j;  1  =  1»  2»  •••  *5*  i  =  2»  •••  »mK  and 

j  =  1,  2,  ,m) .  This  strategy  for  monitoring  production  is  analogous  to  the 

cumulative  result  criterion  of  Cone  and  Dodge  (1964)  in  lot  inspection,  where 
sampling  is  terminated  when  the  cumulative  result  from  a  starting  point  exceeds 
its  criterion. 
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«9.  KEY  WORDS  (Continue  on  reverie  aide  it  nvee  emery  end  Identity  b\ ■  block  number) 

Monitoring  production  processes,  variable  sampling  rates,  Markov  chains, 
number  of  inspections,  number  of  adjustments  to  production. 


JO.  AOSTRACT  (Continue  on  revoreo  aide  If  neceeesry  end  Identity  by  block  number) 

Sampling  procedures  are  considered  for  monitoring  the  output  of  a  pro¬ 
duction  process.  Sampling  rates  are  allowed  to  depend  on  one  of  several 
levels  of  acceptance  or  rejection,  sampling  less  frequently  when  the  process 
is  in  control.  The  problem  is  formulated  as  a  simple  Markov  process  whose 
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properties  yield  the  expected  values  and  the  variances  of  the  number  of  samples 
and  the  numbers  of  various  types  of  adjustments  to  the  process.  Computations 
are  given  in  support  of  the  economic  design  of  variable  sampling  policies  of 
this' type. 
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